EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP-2013-004 
LHCb-PAPER-20 12-045 
January 24, 2013 



Analysis of the resonant components 



Interpretation of CP violation measurements using charmonium decays, in both 
the B° and B® systems, can be subject to changes due to "penguin" type diagrams. 
These effects can be investigated using measurements of the Cabibbo-suppressed 
B° — > J/ipTT + ir~ decays. The final state composition of this channel is investigated 
using a l.Ofb -1 sample of data produced in 7 TeV pp collisions at the LHC and 
collected by the LHCb experiment. A modified Dalitz plot analysis is performed using 
both the invariant mass spectra and the decay angular distributions. An improved 
measurement of the B° — > J/ipir + ir~ branching fraction of (3.97±0.09±0.11±0.16) x 
10 -5 is reported where the first uncertainty is statistical, the second is systematic 
and the third is due to the uncertainty of the branching fraction of the decay B — > 
J/ipK~ used as a normalization channel. Significant production of /o(500) and p(770) 
resonances is found in the substructure of the J/tp-K + ir~ final state, and this indicates 
that they are viable final states for CP violation studies. In contrast evidence for the 
/o(980) resonance is not found. This allows us to establish the first upper limit on the 
branching fraction product B (B° J/tp f (980)) xB (/o(980) -> vr+vr-) < l.lxlfr 6 , 
leading to an upper limit on the absolute value of the mixing angle of the /o(980) 
with the /o(500) of less than 31°, both at 90% confidence level. 



m 




The LHCb collaboration!!] 



Abstract 



Submitted to Physical Review D 



© CERN on behalf of the LHCb collaboration, license CC-BY-3.0 



^Authors are listed on the following pages. 



LHCb collaboration 



R. Aaij 38 , C. Abellan Beteta 33 -™, A. Adametz 11 , B. Adeva 34 , M. Adinolfi 43 , C. Adrover 6 , 

A. Affolder 49 , Z. Ajaltouni 5 , J. Albrecht 9 , F. Alessio 35 , M. Alexander 48 , S. Ali 38 , G. Alkhazov 27 , 

P. Alvarez Cartelle 34 , A. A. Alves Jr 22 ' 35 , S. Amato 2 , Y. Amhis 7 , L. Anderlini 17 '^, J. Anderson 37 , 

R. Andreassen 57 , R.B. Appleby 51 , O. Aquines Gutierrez 10 , F. Archilli 18 , A. Artamonov 32 , 

M. Artuso 53 , E. Aslanides 6 , G. Auriemma 22 '" 1 , S. Bachmann 11 , J.J. Back 45 , C. Baesso 54 , 

V. Balagura 28 , W. Baldini 16 , R.J. Barlow 51 , C. Barschel 35 , S. Barsuk 7 , W. Barter 44 , 

Th. Bauer 38 , A. Bay 36 , J. Beddow 48 , I. Bediaga 1 , S. Belogurov 28 , K. Belous 32 , I. Belyaev 28 , 

E. Ben-Haim 8 , M. Benayoun 8 , G. Bencivenni 18 , S. Benson 47 , J. Benton 43 , A. Berezhnoy 29 , 

R. Bernet 37 , M.-O. Bettler 44 , M. van Beuzekom 38 , A. Bien 11 , S. Bifani 12 , T. Bird 51 , 

A. Bizzeti 17 '' 1 , P.M. Bj0rnstad 51 , T. Blake 35 , F. Blanc 36 , C. Blanks 50 , J. Blouw 11 , S. Blusk 53 , 

A. Bobrov 31 , V. Bocci 22 , A. Bondar 31 , N. Bondar 27 , W. Bonivento 15 , S. Borghi 51 , A. Borgia 53 , 

T.J.V. Bowcock 49 , E. Bowen 37 , C. Bozzi 16 , T. Brambach 9 , J. van den Brand 39 , J. Bressieux 36 , 

D. Brett 51 , M. Britsch 10 , T. Britton 53 , N.H. Brook 43 , H. Brown 49 , I. Burducea 26 , A. Bursche 37 , 

J. Buytaert 35 , S. Cadeddu 15 , O. Callot 7 , M. Calvi 20 ^', M. Calvo Gomez 33 '™, A. Camboni 33 , 

P. Campana 18 - 35 , A. Carbone 14 ' c , G. Carboni 21 ' fc , R. Cardinale 19 '\ A. Cardini 15 , 

H. Carranza-Mejia 47 , L. Carson 50 , K. Carvalho Akiba 2 , G. Casse 49 , M. Cattaneo 35 , Ch. Cauet 9 , 

M. Charles 52 , Ph. Charpentier 35 , P. Chen 3 ' 36 , N. Chiapolini 37 , M. Chrzaszcz 23 , K. Ciba 35 , 

X. Cid Vidal 34 , G. Ciezarek 50 , P.E.L. Clarke 47 , M. Clemencic 35 , H.V. Cliff 44 , J. Closier 35 , 

C. Coca 26 , V. Coco 38 , J. Cogan 6 , E. Cogneras 5 , P. Collins 35 , A. Comerma-Montells 33 , 

A. Contu 15 , A. Cook 43 , M. Coombes 43 , G. Corti 35 , B. Couturier 35 , G.A. Cowan 36 , D. Craik 45 , 
S. Cunliffe 50 , R. Currie 47 , C. D'Ambrosio 35 , P. David 8 , P.N.Y. David 38 , I. De Bonis 4 , 
K. De Bruyn 38 , S. De Capua 51 , M. De Cian 37 , J.M. De Miranda 1 , L. De Paula 2 , W. De Silva 57 , 
P. De Simone 18 , D. Decamp 4 , M. Deckenhoff 9 , H. Degaudenzi 36 ' 35 , L. Del Buono 8 , C. Deplano 15 , 

D. Derkach 14 , O. Deschamps 5 , F. Dettori 39 , A. Di Canto 11 , J. Dickens 44 , H. Dijkstra 35 , 
P. Diniz Batista 1 , M. Dogaru 26 , F. Domingo Bonal 33 '™, S. Donleavy 49 , F. Dordei 11 , 

A. Dosil Suarez 34 , D. Dossett 45 , A. Dovbnya 40 , F. Dupertuis 36 , R. Dzhelyadin 32 , A. Dziurda 23 , 
A. Dzyuba 27 , S. Easo 46 ' 35 , U. Egede 50 , V. Egorychev 28 , S. Eidelman 31 , D. van Eijk 38 , 
S. Eisenhardt 47 , U. Eitschberger 9 , R. Ekelhof 9 , L. Eklund 48 , I. El Rifai 5 , Ch. Elsasser 37 , 
D. Elsby 42 , A. Falabella 14 ' 6 , C. Farber 11 , G. Fardell 47 , C. Farinelli 38 , S. Farry 12 , V. Fave 36 , 

D. Ferguson 47 , V. Fernandez Albor 34 , F. Ferreira Rodrigues 1 , M. Ferro-Luzzi 35 , S. Filippov 30 , 
C. Fitzpatrick 35 , M. Fontana 10 , F. Fontanelli 19 '*, R. Forty 35 , O. Francisco 2 , M. Frank 35 , 

C. Frei 35 , M. Frosini 17 '^, S. Furcas 20 , E. Furfaro 21 , A. Gallas Torreira 34 , D. Galli 14 ' c , 

M. Gandelman 2 , P. Gandini 52 , Y. Gao 3 , J. Garofoli 53 , P. Garosi 51 , J. Garra Tico 44 , 

L. Garrido 33 , C. Gaspar 35 , R. Gauld 52 , E. Gersabeck 11 , M. Gersabeck 51 , T. Gershon 45 ' 35 , 

Ph. Ghez 4 , V. Gibson 44 , V.V. Gligorov 35 , C. Gobel 54 , D. Golubkov 28 , A. Golutvin 50 ' 28 ' 35 , 

A. Gomes 2 , H. Gordon 52 , M. Grabalosa Gandara 5 , R. Graciani Diaz 33 , L.A. Granado Cardoso 35 , 

E. Grauges 33 , G. Graziani 17 , A. Grecu 26 , E. Greening 52 , S. Gregson 44 , O. Griinberg 55 , B. Gui 53 , 

E. Gushchin 30 , Yu. Guz 32 , T. Gys 35 , C. Hadjivasiliou 53 , G. Haefeli 36 , C. Haen 35 , S.C. Haines 44 , 
S. Hall 50 , T. Hampson 43 , S. Hansmann-Menzemer 11 , N. Harnew 52 , S.T. Harnew 43 , J. Harrison 51 , 
P.F. Harrison 45 , T. Hartmann 55 , J. He 7 , V. Heijne 38 , K. Hennessy 49 , P. Henrard 5 , 

J.A. Hernando Morata 34 , E. van Herwijnen 35 , E. Hicks 49 , D. Hill 52 , M. Hoballah 5 , 

C. Hombach 51 , P. Hopchev 4 , W. Hulsbergen 38 , P. Hunt 52 , T. Huse 49 , N. Hussain 52 , 

D. Hutchcroft 49 , D. Hynds 48 , V. Iakovenko 41 , P. Ilten 12 , R. Jacobsson 35 , A. Jaeger 11 , E. Jans 38 , 

F. Jansen 38 , P. Jaton 36 , F. Jing 3 , M. John 52 , D. Johnson 52 , C.R. Jones 44 , B. Jost 35 , M. Kaballo 9 , 



ii 



S. Kandybei 40 , M. Karacson 35 , T.M. Karbach 35 , I.R. Kenyon 42 , U. Kerzel 35 , T. Ketel 39 , 

A. Keune 36 , B. Khanji 20 , O. Kochebina 7 , I. Komarov 36 ' 29 , R.F. Koopman 39 , P. Koppenburg 38 , 

M. Korolev 29 , A. Kozlinskiy 38 , L. Kravchuk 30 , K. Kreplin 11 , M. Kreps 45 , G. Krocker 11 , 

P. Krokovny 31 , F. Kruse 9 , M. Kucharczyk 20 ' 23 '- 5 ; , V. Kudryavtsev 31 , T. Kvaratskheliya 28 ' 35 , 

V.N. La Thi 36 , D. Lacarrere 35 , G. Lafferty 51 , A. Lai 15 , D. Lambert 47 , R.W. Lambert 39 , 

E. Lanciotti 35 , G. Lanfranchi 18,35 , C. Langenbruch 35 , T. Latham 45 , C. Lazzeroni 42 , R. Le Gac 6 , 

J. van Leerdam 38 , J.-P. Lees 4 , R. Lefevre 5 , A. Leflat 29 ' 35 , J. Lefrangois 7 , O. Leroy 6 , Y. Li 3 , 

L. Li Gioi 5 , M. Liles 49 , R. Lindner 35 , C. Linn 11 , B. Liu 3 , G. Liu 35 , J. von Loeben 20 , J.H. Lopes 2 , 

E. Lopez Asamar 33 , N. Lopez-March 36 , H. Lu 3 , J. Luisier 36 , H. Luo 47 , F. Machefert 7 , 

I.V. Machikhiliyan 4 ' 28 , F. Maciuc 26 , O. Maev 27 ' 35 , S. Malde 52 , G. Manca 15 ' d , G. Mancinelli 6 , 

N. Mangiafave 44 , U. Marconi 14 , R. Marki 36 , J. Marks 11 , G. Martellotti 22 , A. Martens 8 , 

L. Martin 52 , A. Martin Sanchez 7 , M. Martinelli 38 , D. Martinez Santos 39 , D. Martins Tostes 2 , 

A. Massafferri 1 , R. Matev 35 , Z. Mathe 35 , C. Matteuzzi 20 , M. Matveev 27 , E. Maurice 6 , 

A. Mazurov 16 ' 30 ' 35 ' 6 , J. McCarthy 42 , R. McNulty 12 , B. Meadows 57,52 , F. Meier 9 , M. Meissner 11 , 
M. Merk 38 , D.A. Milanes 8 , M.-N. Minard 4 , J. Molina Rodriguez 54 , S. Monteil 5 , D. Moran 51 , 

P. Morawski 23 , R. Mountain 53 , I. Mous 38 , F. Muheim 47 , K. Miiller 37 , R. Muresan 26 , B. Muryn 24 , 

B. Muster 36 , P. Naik 43 , T. Nakada 36 , R. Nandakumar 46 , I. Nasteva 1 , M. Needham 47 , 
N. Neufeld 35 , A.D. Nguyen 36 , T.D. Nguyen 36 , C. Nguyen-Mau 36 '°, M. Nicol 7 , V. Mess 5 , 
R. Met 9 , N. Nikitin 29 , T. Nikodem 11 , S. Nisar 56 , A. Nomerotski 52 , A. Novoselov 32 , 

A. Oblakowska-Mucha 24 , V. Obraztsov 32 , S. Og gero 38 , S. Ogilvy 48 , O. Okhrimenko 41 , 

R. 01deman 15 ' d ' 35 , M. Orlandea 26 , J.M. Otalora Goicochea 2 , P. Owen 50 , B.K. Pal 53 , 

A. Palano 13 ' 6 , M. Palutan 18 , J. Panman 35 , A. Papanestis 46 , M. Pappagallo 48 , C. Parkes 51 , 

C. J. Parkinson 50 , G. Passaleva 17 , G.D. Patel 49 , M. Patel 50 , G.N. Patrick 46 , C. Patrignani 19 '\ 
C. Pavel-Nicorescu 26 , A. Pazos Alvarez 34 , A. Pellegrino 38 , G. Penso 22 '', M. Pepe Altarelli 35 , 
S. Perazzini 14,c , D.L. Perego 20 '- 5 , E. Perez Trigo 34 , A. Perez-Calero Yzquierdo 33 , P. Perret 5 , 
M. Perrin-Terrin 6 , G. Pessina 20 , K. Petridis 50 , A. Petrolini 19 '*, A. Phan 53 , 

E. Picatoste Olloqui 33 , B. Pietrzyk 4 , T. Pilar 45 , D. Pinci 22 , S. Playfer 47 , M. Plo Casasus 34 , 

F. Polci 8 , G. Polok 23 , A. Poluektov 45 ' 31 , E. Polycarpo 2 , D. Popov 10 , B. Popovici 26 , 

C. Potterat 33 , A. Powell 52 , J. Prisciandaro 36 , V. Pugatch 41 , A. Puig Navarro 36 , W. Qian 4 , 
J.H. Rademacker 43 , B. Rakotomiaramanana 36 , M.S. Rangel 2 , I. Raniuk 40 , N. Rauschmayr 35 , 

G. Raven 39 , S. Redford 52 , M.M. Reid 45 , A.C. dos Reis 1 , S. Ricciardi 46 , A. Richards 50 , 
K. Rinnert 49 , V. Rives Molina 33 , D.A. Roa Romero 5 , P. Robbe 7 , E. Rodrigues 51 , 

P. Rodriguez Perez 34 , G.J. Rogers 44 , S. Roiser 35 , V. Romanovsky 32 , A. Romero Vidal 34 , 

J. Rouvinet 36 , T. Ruf 35 , H. Ruiz 33 , G. Sabatino 22 < fc , J.J. Saborido Silva 34 , N. Sagidova 27 , 

P. Sail 48 , B. Saitta 15 ' d , C. Salzmann 37 , B. Sanmartin Sedes 34 , M. Sannino 19 '*, R. Santacesaria 22 , 

C. Santamarina Rios 34 , E. Santovetti 21 ' fc , M. Sapunov 6 , A. Sarti 18 '', C. Satriano 22 '™, A. Satta 21 , 

M. Savrie 16 ' 6 , D. Savrina 28 ' 29 , P. Schaack 50 , M. Schiller 39 , H. Schindler 35 , S. Schleich 9 , 

M. Schlupp 9 , M. Schmelling 10 , B. Schmidt 35 , O. Schneider 36 , A. Schopper 35 , M.-H. Schune 7 , 

R. Schwemmer 35 , B. Sciascia 18 , A. Sciubba 18 ' Z , M. Seco 34 , A. Semennikov 28 , K. Senderowska 24 , 

I. Sepp 50 , N. Serra 37 , J. Serrano 6 , P. Seyfert 11 , M. Shapkin 32 , I. Shapoval 40 ' 35 , P. Shatalov 28 , 

Y. Shcheglov 27 , T. Shears 49 ' 35 , L. Shekhtman 31 , O. Shevchenko 40 , V. Shevchenko 28 , A. Shires 50 , 

R. Silva Coutinho 45 , T. Skwarnicki 53 , N.A. Smith 49 , E. Smith 52 ' 46 , M. Smith 51 , K. Sobczak 5 , 

M.D. Sokoloff 57 , F.J.P. Soler 48 , F. Soomro 18 ' 35 , D. Souza 43 , B. Souza De Paula 2 , B. Spaan 9 , 

A. Sparkes 47 , P. Spradlin 48 , F. Stagni 35 , S. Stahl 11 , O. Steinkamp 37 , S. Stoica 26 , S. Stone 53 , 

B. Storaci 37 , M. Straticiuc 26 , U. Straumann 37 , V.K. Subbiah 35 , S. Swientek 9 , V. Syropoulos 39 , 
M. Szczekowski 25 , P. Szczypka 36 ' 35 , T. Szumlak 24 , S. T'Jampens 4 , M. Teklishyn 7 , 



iii 



E. Teodorescu 26 , F. Teubert 35 , C. Thomas 52 , E. Thomas 35 , J. van Tilburg 11 , V. Tisserand 4 , 

M. Tobin 37 , S. Tolk 39 , D. Tonelli 35 , S. Topp-Joergensen 52 , N. Torr 52 , E. Tournefier 4 - 50 , 

S. Tourneur 36 , M.T. Tran 36 , M. Tresch 37 , A. Tsaregorodtsev 6 , P. Tsopelas 38 , N. Tuning 38 , 

M. Ubeda Garcia 35 , A. Ukleja 25 , D. Urner 51 , U. Uwer 11 , V. Vagnoni 14 , G. Valenti 14 , 

R. Vazquez Gomez 33 , P. Vazquez Regueiro 34 , S. Vecchi 16 , J.J. Velthuis 43 , M. Veltri 17,9 , 

G. Veneziano 36 , M. Vesterinen 35 , B. Viaud 7 , D. Vieira 2 , X. Vilasis-Cardona 33 '™, A. Vollhardt 37 , 

D. Volyanskyy 10 , D. Voong 43 , A. Vorobyev 27 , V. Vorobyev 31 , C. Vofi 55 , H. Voss 10 , R. Waldi 55 , 

R. Wallace 12 , S. Wandernoth 11 , J. Wang 53 , D.R. Ward 44 , N.K. Watson 42 , A.D. Webber 51 , 

D. Websdale 50 , M. Whitehead 45 , J. Wicht 35 , J. Wiechczynski 23 , D. Wiedner 11 , L. Wiggers 38 , 

G. Wilkinson 52 , M.P. Williams 45 ' 46 , M. Williams 50 ^, F.F. Wilson 46 , J. Wishahi 9 , M. Witek 23 , 

S.A. Wotton 44 , S. Wright 44 , S. Wu 3 , K. Wyllie 35 , Y. Xie 47 - 35 , F. Xing 52 , Z. Xing 53 , Z. Yang 3 , 

R. Young 47 , X. Yuan 3 , O. Yushchenko 32 , M. Zangoli 14 , M. Zavertyaev 10,a , F. Zhang 3 , 

L. Zhang 53 , W.C. Zhang 12 , Y. Zhang 3 , A. Zhelezov 11 , A. Zhokhov 28 , L. Zhong 3 , A. Zvyagin 35 . 

1 Centro Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 

2 Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 

3 Center for High Energy Physics, Tsinghua University, Beijing, China 
4 LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

5 Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
e CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
7 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

S LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 
9 Fakultdt Physik, Technische Universitdt Dortmund, Dortmund, Germany 
10 'Max- Planck- Institut fur Kernphysik (MPIK), Heidelberg, Germany 

11 Physikalisches Institut, Ruprecht-Karls- Universitdt Heidelberg, Heidelberg, Germany 

12 School of Physics, University College Dublin, Dublin, Ireland 



13 



Sezione INFN di 



Bari, Bari, Italy 



14 Sezione INFN di Bologna, Bologna, Italy 

15 Sezione INFN di Cagliari, Cagliari, Italy 

16 Sezione INFN di Ferrara, Ferrara, Italy 

17 Sezione INFN di Firenze, Firenze, Italy 

18 Laboratori Nazionali dell 'INFN di Frascati, Frascati, Italy 

19 Sezione INFN di Genova, Genova, Italy 

20 'Sezione INFN di Milano Bicocca, Milano, Italy 

21 Sezione INFN di Roma Tor Vergata, Roma, Italy 

22 Sezione INFN di Roma La Sapienza, Roma, Italy 

23 Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland 
24 AGH University of Science and Technology, Krakow, Poland 
25 National Center for Nuclear Research (NCBJ), Warsaw, Poland 

26 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

27 Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

28 Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

29 Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 

30 Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 

31 Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 

32 Institute for High Energy Physics (IHEP), Protvino, Russia 

33 Universitat de Barcelona, Barcelona, Spain 

34 Universidad de Santiago de Compostela, Santiago de Compostela, Spain 

35 European Organization for Nuclear Research (CERN), Geneva, Switzerland 
36 Ecole Polytechnique Federate de Lausanne (EPFL), Lausanne, Switzerland 



iv 



37 Physik-Institut, Universitat Zurich, Zurich, Switzerland 

38 Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 

39 Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, The 

Netherlands 

40 NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 

41 Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 

42 University of Birmingham, Birmingham, United Kingdom 

43 H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

44 Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 

45 Department of Physics, University of Warwick, Coventry, United Kingdom 
46 STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

47 School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
48 School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 
49 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 
50 Imperial College London, London, United Kingdom 

51 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 

52 Department of Physics, University of Oxford, Oxford, United Kingdom 

53 Syracuse University, Syracuse, NY, United States 

54 Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to 2 

55 Institut fur Physik, Universitat Rostock, Rostock, Germany, associated to 11 

56 Institute of Information Technology, COMSATS, Lahore, Pakistan, associated to 53 

57 University of Cincinnati, Cincinnati, OH, United States, associated to 53 

a P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 



b Universita 
c Universita di 
d Universita di 
e Universita di 



Bari, Bari, Italy 
Bologna, Bologna, Italy 
Cagliari, Cagliari, Italy 
Ferrara, Ferrara, Italy 



f Universita di Firenze, Firenze, Italy 
9 Universita di Urbino, Urbino, Italy 

Universita di Modena e Reggio Emilia, Modena, Italy 
1 Universita di Genova, Genova, Italy 
J Universita di Milano Bicocca, Milano, Italy 
k Universita di Roma Tor Vergata, Roma, Italy 
1 Universita di Roma La Sapienza, Roma, Italy 
m Universita delta Basilicata, Potenza, Italy 

n LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain 
"Hanoi University of Science, Hanoi, Viet Nam 

p Massachusetts Institute of Technology, Cambridge, MA, United States 



V 



1 Introduction 



CP violation measurements using neutral B meson decays into J/0 mesons are of prime 
importance both for determinations of Standard Model (SM) parameters and searching 
for physics beyond the SM. In the case of B° decays, the final state J/tpK® is the most 
important for measuring sin 2/3 [1], while in the case of B° s decays, used to measure <f> s , 
only the final states J/00 (2]-[4j, and J/07r + 7r~ [5] have been used so far, where the largest 
component of the latter is J/0 /o (980) pi. The decay rate for these J/ip modes is dominated 
by the color-suppressed tree level diagram, an example of which is shown for B° decays in 
Fig. [IJa), while penguin processes, an example of which is shown in Fig. [IJb), are expected 
to be suppressed. Theoretical predictions on the effects of such "penguin pollution" vary 
widely for both B° and B° s decays [7j, so it is incumbent upon experimentalists to limit 
possible changes in the value of the CP violating angles measured using other decay modes. 




Figure 1: (a) Tree level and (b) penguin diagram examples for B° decays into J/ipK®. 



The decay B° — > J/ip7t + TT~ can occur via a Cabibbo suppressed tree level diagram, 
shown in Fig. [2^a), or via several penguin diagrams. An example is shown in Fig. [2^b), 
while others are illustrated in Ref. [81. These decays are interesting because they can also 
be used to measure or limit the amount of penguin pollution. The advantage in using 
the decay B° — > J/?/>7r + 7r~ arises because the relative amount of pollution is larger. In 
the allowed decays, e.g. J? — > J/ipK®, the penguin amplitude is multiplied by a factor of 
\ 2 Re 1 ^, where A is the sine of the Cabibbo angle (~ 0.22), while in the suppressed decays 
the factor becomes R'e^ , where R and R', and and 0' are expected to be similar in 
size J8j. A similar study uses the decay B® — > J/ipK® [§]. 

CP violation measurements in the J/if)7r + 7T~ mode utilizing B° — B° mixing determine 
sin2/3 eff which can be compared to the well measured sin 2/3. Differences can be used to 
estimate the magnitude of penguin effects. Knowledge of the final state structure is the 
first step in this program. Such measurements on sin 2/3 eff have been attempted in the B° 




Figure 2: (a) Tree level and (b) penguin diagram for B° decays into J/ifj7i + 7i 
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system by using the J/tpir final state 10 



In order to ascertain the viability of such CP violation measurements we perform a full 
"Dalitz like" analysis of the final state. Regions in 7r + 7r~ mass that correspond to spin-0 
final states would be CP eigenstates. Final states containing vector resonances, such as the 
p(770) can be analyzed in a similar manner as was done for the decay B° s — > J/ipcj) j^jjl]. 

It is also of interest to search for the /o(980) contribution and to obtain information 
concerning the mixing angle between the /o(980) and the /o(500)Q partners in the scalar 
nonet, as the latter should couple strongly to the dd system. Branching fractions for 
B° — > J/il)n + 'K~ and J/tp p° have previously been measured by the BaBar collaboration 11 



In this paper the J/vpic + and tt + -k~ mass spectra and decay angular distributions are 
used to determine the resonant and non- resonant components. This differs from a classical 
Dalitz plot analysis 12 because one of the particles in the final state, the J/ip meson, has 
spin-1 and its three decay amplitudes must be considered. We first show that there are 
no evident structures in the <//?/> 7r + invariant mass, and then model the ir + n~ invariant 
mass with a series of resonant and non- resonant amplitudes. The data are then fitted with 
the coherent sum of these amplitudes. We report on the resonant structure and the CP 
content of the final state. 



2 Data sample and selection requirements 

The data sample consists of 1.0 fb _1 of integrated luminosity collected with the LHCb 



detector 13 using pp collisions at a center-of-mass energy of 7 TeV. The detector is a single- 
arm forward spectrometer covering the pseudorapidity range 2 < rj < 5, designed for the 
study of particles containing b or c quarks. Components include a high precision tracking 
system consisting of a silicon-strip vertex detector surrounding the pp interaction region, 
a large-area silicon-strip detector located upstream of a dipole magnet with a bending 
power of about 4 Tm, and three stations of silicon-strip detectors and straw drift-tubes 
placed downstream. The combined tracking system has a momentum^] resolution Ap/p 
that varies from 0.4% at 5 GeV to 0.6% at 100 GeV, and an impact parameter resolution 
of 20 \mi for tracks with large transverse momentum (jot) with respect to the proton 
beam direction. Charged hadrons are identified using two ring-imaging Cherenkov (RICH) 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a system composed of alternating layers of 
iron and multiwire proportional chambers. The trigger consists of a hardware stage, based 
on information from the calorimeter and muon systems, followed by a software stage that 
applies a full event reconstruction |14|. 

Events are triggered by a J/ip — > decay, requiring two identified muons with 

opposite charge, pri^) greater than 500 MeV, an invariant mass within 120 MeV of 
the J/ip mass [15], and form a vertex with a fit \ 2 l ess than 16. After applying these 



1 This particle has been identified previously as the Jo (600) or a resonance. 
2 We work in units where c = 1. 
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requirements, there is a large J/ip signal over a small background 16 . Only candidates 



with dimuon invariant mass between — 48MeV and +43 MeV relative to the observed 
J/ip mass peak are selected, corresponding a window of about ±3cr. The requirement is 
asymmetric because of final state electromagnetic radiation. The two muons subsequently 
are kinematically constrained to the known J/ijj mass. 

Other requirements are imposed to isolate B° candidates with high signal yield and 
minimum background. This is accomplished by combining the J/ip — > /i + pr candidate 
with a pair of pion candidates of opposite charge, and then testing if all four tracks form a 
common decay vertex. Pion candidates are each required to have pt greater than 250 MeV, 
and the scalar sum of the two transverse momenta, Pt{^ + ) + Pt{tt~), must be larger than 
900 MeV. The impact parameter (IP) is the distance of closest approach of a track to the 
primary vertex (PV). To test for inconsistency with production at the PV, the IP y 2 is 
computed as the difference between the y 2 of the PV reconstructed with and without the 
considered track. Each pion must have an IP y 2 greater than 9. Both pions must also 
come from a common vertex with an acceptable y 2 an d form a vertex with the J/ip with a 
y 2 per number of degrees of freedom (ndf) less than 10 (here ndf equals five). Pion and 
kaon candidates are positively identified using the RICH system. Cherenkov photons are 
matched to tracks, the emission angles of the photons compared with those expected if 
the particle is an electron, pion, kaon or proton, and a likelihood is then computed. The 
particle identification makes use of the logarithm of the likelihood ratio comparing two 
particle hypotheses (DLL). For pion selection we require DLL(7r — K) > —10. 

The four-track B° candidate must have a flight distance of more than 1.5 mm, where the 
average decay length resolution is 0.17 mm. The angle between the combined momentum 
vector of the decay products and the vector formed from the positions of the PV and the 
decay vertex (pointing angle) is required to be less than 2.5°. 

Events satisfying this preselection are then further filtered using a multivariate analyzer 



based on a Boosted Decision Tree (BDT) technique 17 . The BDT uses six variable 



that are chosen in a manner that does not introduce an asymmetry between either the 
two muons or the two pions. They are the minimum DLL(/z — 7r) of the fi + and [T , the 
minimum p^ of the 7r + and 7r~, the minimum of the IP y 2 of the 7r + and 7r~, the B° 
vertex y 2 , the B° pointing angle, and the B° flight distance. There is discrimination power 
between signal and background in all of these variables, especially the B° vertex y 2 . 

The background sample used to train the BDT consists of the events in the B° mass 
sideband having 5566 < m( <//?/> 7r + 7r~) < 5616 MeV. The signal sample consists of two 
million B° — > J/i[>{— > /U + /i~)7r + 7r~ Monte Carlo simulated events that are generated 



uniformly in phase space, using Pythia |18| with a special LHCb parameter tune 19 



and the LHCb detector simulation based on Geant4 20 described in Ref. 21 . Separate 
samples are used to train and test the BDT. The distributions of the BDT classifier for 
signal and background are shown in Fig. [3j To minimize a possible bias on the signal 
acceptance due to the BDT, we choose a relatively loose requirement of the BDT classifier 
> 0.05 which has a 96% signal efficiency and a 92% background rejection rate. 

The invariant mass of the selected J/ip7i + ii^ combinations, where the dimuon pair is 
constrained to have the J/ip mass, is shown in Fig. g There are signal peaks at both 
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Figure 3: Distributions of the BDT classifier for both training and test samples of J/ip-K + 7T~ 
signal and background events. The signal samples are from simulation and the background 
samples are from data. 
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Figure 4: Invariant mass of J/tfjir + 7T~ combinations. The data are fitted with a double- 
Gaussian signal and several background functions. The (red) solid double-Gaussian function 
centered at 5280 MeV is the B° signal, the (brown) dotted line shows the combinatorial 
background, the (green) short-dashed shows the B~ background, the (purple) dot-dashed 
line shows the contribution of B° s — > J/ip7i + ii~ decays, the (black) dot-long dashed is the 
sum of B° s — > J/if)r/(—> pj) and B° s — > J/i/;(f>(—> 7r + 7r~7r°) backgrounds, the (light blue) 
long-dashed is the B° — > J/ifjK~ir + reflection, and the (blue) solid line is the total. 
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the B° s and B° masses on top of the background. Double-Gaussian functions are used 
to fit both signal peaks. They differ only in their mean values, which are determined by 
the data. The core Gaussian width is also allowed to vary, while the fraction and width 
ratio of the second Gaussian is fixed to that obtained in the fit of B° s — > J/ip(p events. 
(The details of the fit are given in Ref. [6].) Other components in the fit model take 
into account background contributions. One source is from B~ — > J/ip K~ decays, which 
contributes when the K~ is misidentifed as a tt~ and then combined with a random 7r + ; 
the smaller J/ipTt~ mode contributes when it is combined with a random 7r + . The next 
source contains B° s — > Jjipvj '(— )■ p'j) and B° s — > J/ip(p(— > 7r + 7r~7r°) decays where the 7 and 
the 7T° are ignored respectively. Finally there is a B° — > J/ipK~ir + reflection where the K~ 
is misidentified as ir~ . Here and elsewhere charged conjugated modes are included when 
appropriate. The exponential combinatorial background shape is taken from same-sign 
combinations, that are the sum of J/ipir + 7i + and J/ip 7r~7r~ candidates. The shapes of the 
other components are taken from the simulation with their normalizations allowed to vary. 
The fit gives 5287 ± 112 signal and 3212 ± 80 background candidates within ±20 MeV of 
the B° mass peak, where a K® veto, discussed later, is applied. 

We use the well measured B~ —t J/ip K~ mode as a normalization channel to determine 
the branching fractions. To minimize the systematic uncertainty from the BDT selection, 
we employ a similar selection on B~ — > J/ipK~ decays after requiring the same pre- 
selection except for particle identification criteria on the K~ candidates. Similar variables 
are used for the BDT except that the variables describing the combination of tt + and 
7T~ in the J/ip7r + 7r~ final state are replaced by ones describing the K~ meson. For BDT 
training, the signal sample uses simulated events and the background sample consists of 
the data events in the region 5400 < m(J/ipK~) < 5450 MeV. The resulting invariant 
mass distribution of the candidates satisfying BDT classifier > 0.05 is shown in Fig. [5j 
Fitting the distribution with a double-Gaussian function for the signal and linear function 
for the background gives 350,727 ± 633 signal and 4756 ± 103 background candidates 
within ±20 MeV of the B~ mass peak. 



3 Analysis formalism 



We apply a formalism similar to that used in Belle's analysis 22 of B —t K~ir + x c i 
decays and later used in LHCb's analysis of B° s — > J/ipir + -ir~ decays [6|. The decay 
B° — > J/tpn + 7i~ , with J/ip — > can be described by four variables. These are taken 

to be the invariant mass squared of J/ipTi + (s\2 = m 2 (J/^7r + )), the invariant mass squared 
of 7r + 7r~ (S23 = in 2 (n + 7i~)) , where we use label 1 for J/ip , 2 for 7r + and 3 for 7r~, the 
J/ip helicity angle {Oj/^,), which is the angle of the fi + in the J/ip rest frame with respect 
to the J/ip direction in the B° rest frame, and the angle between the J/ip and 7r + 7r~ 
decay planes (x) m the B° rest frame. To improve the resolution of these variables we 
perform a kinematic fit constraining the B° and J/ip masses to their nominal values 15 



and recompute the final state momenta. To simplify the probability density function, we 
analyze the decay process after integrating over Xi which eliminates several interference 
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Figure 5: Invariant mass of J/%pK~ combinations. The data points are fitted with a 
double-Gaussian function for signal and a linear function for background. The dotted line 
shows the background, and the (blue) solid line is the total. 



terms. 



3.1 The decay model for B° — > J/i/?7r + 7r 

The overall probability density function (PDF) given by the sum of signal, S, and 
background, B, functions is 

F(Si2,S 2 3,0j/i>) = -T^g(-Sl2, S 2 3,0j/ip)S(si2, S 2 3,0j/^) + ~ 7 r B(s 12 , S 23 , Ojfy ) , (1) 

yVgig vVbkg 

where / s j g is the fraction of the signal in the fitted region and e is the detection efficiency. 
The fraction of the signal is obtained from the mass fit and is fixed for the subsequent 
analysis. The normalization factors are given by 

Mag = J £(s 12 ,s 23 ,9j/^)S(si 2 , s 23 ,9 J/i(> ) ds 12 ds 23 d cos 9 j/^, 

TVbkg = J B(s 12 ,s 23 , 9 jfy) ds 12 ds 23 d cos 9 j^. (2) 



The event distribution for m 2 (n + n ) versus m 2 (J/ip7i + ) in Fig. [6] shows obvious 
structure in m 2 (7r + 7r~). To investigate if there are visible exotic structures in the J/i/j7t + 



system as claimed in similar decays 23 , we examine the J/ip ir + mass distribution shown in 
Fig. (7b a). No resonant effects are evident. Figure [7] (b) shows the tt + -k~ mass distribution. 
There is a clear peak at the p(770) region, a small bump around 1250 MeV, but no evidence 
for the /q(980) resonance. The favored B° — > J/tpK® decay is mostly rejected by the B° 
vertex x 2 selection, but about 150 such events remain. We eliminate them by excluding 
the candidates that have |m(7r + 7r~) — m K o \ < 25 MeV, where m K o is the K® mass |15j . 
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Figure 6: Distribution of m 2 (ir + 7i ) versus m 2 ( J/ip tt + ) for 5° candidate decays within 
±20 MeV of the B° mass. 
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Figure 7: Distribution of (a) m(J/ip7T + ) and (b) m(7r + 7r~) for S° — > J/ip7T + 7T~ candidate 
decays within ±20 MeV of B° mass shown with the solid line. The (red) points with error 
bars show the background contribution determined from m( J/ip tt + 7t~) fits performed in 
each bin. 
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3.1.1 The signal function 

The signal function for B° is taken to be the coherent sum over resonant states that can 
decay into 7r + 7r+ plus a possible non-resonant S-wave contribution^] 



S(si2,S23,0jfr) = 



A=0,±1 



Y,^e< i Af{s l2 ,s 2 ^e J ^) 



(3) 



where A^{si2, s 23 , Oj/^p) is the amplitude of the decay via an intermediate resonance Ri 
with helicity A. Each Ri has an asso 
A and a phase Note that the s 
amplitudes for each i are defined as 



with helicity A. Each Ri has an associated amplitude strength a^ 4 for each helicity state 
A and a phase Note that the spin-0 component can only have a A = term. The 



A%(s 12 ,s 23 ,9 J/i} ) = F b Lb) F r Lr) A R (s 23 ) (^-) LB (^=) LR T x Q x (9 J/f ), (4) 



where Pb is the J/ip momentum in the B° rest frame and Pr is the momentum of either 
of the two pions in the dipion rest frame, mg is the E>° mass, F b Lb ^ and F r Lr ^ are the 
B° meson and R resonance Blatt-Weisskopf barrier factors 24 , Lb is the orbital angular 
momentum between the J/tp and 7r + 7r~ system, and Lb is the orbital angular momentum 
in the decay and is equal to the spin of resonance R because pions have spin-0. 

Since the parent B° has spin-0 and the J/tp is a vector, when the tt + tc~ system forms a 
spin-0 resonance, Lb = 1 and Lb = 0. For 7r + 7r~ resonances with non-zero spin, Lb can 
be 0, 1 or 2 (1, 2 or 3) for Lb = 1(2) and so on. We take the lowest Lg as the default and 
consider the other possibilities in the systematic uncertainty. 
The Blatt-Weisskopf barrier factors F b Lb ^ and F r Lr ^ are 

F (0) = h 

F m = v^p, (5) 

Vl + z 



i?(2) 



z z Q +3;zo + 9 



\lz 2 + 3z + 9 



For the B meson z = r 2 P B , where the hadron scale r is taken as 5.0 GeV , and for the R 



resonance z = r 2 P R with r taken as 1.5 GeV 1 25 . In both cases zq = r 2 Pn where Pq is 



the decay daughter momentum calculated at the resonance pole mass. 

The angular term, T x , is obtained using the helicity formalism and is defined as 



T x ^d J X0 (9^), (6) 

where d is the Wigner d-function, J is the resonance spin, 9 nn is the 7r + 7r~ resonance 
helicity angle which is defined as the angle of the tt + in the ir + n~ rest frame with respect 



3 The interference terms between different helicities are zero because we integrate over the angular 
variable x- 
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to the tt + tc direction in the B° rest frame and calculated from the other variables as 

[m 2 (J/ipir + ) -m 2 (J/ipn-)]m(n+TT-) 

cos ^ = ■ (7) 

The J/ip helicity dependent term Q\(0j/tp) is defined as 



®\{0m) = \l sin 2 9ju for helicity = 



\\yj/i>) — V u J/ip 

for helicity = ±1. (8) 



1 + COS 2 dj/^p 



The function Ar{s2^) describes the mass squared shape of the resonance R, that in 
most cases is a Breit-Wigner (BW) amplitude. Complications arise, however, when a new 
decay channel opens close to the resonant mass. The proximity of a second threshold 
distorts the line shape of the amplitude. This happens for the /o(980) resonance because 
the K + K~ decay channel opens. Here we use a Flatte model [26] which is described below. 

The BW amplitude for a resonance decaying into two spin-0 particles, labeled as 2 and 
3, is 

Msn) = — r , x , (9) 

where is the resonance pole mass, T(s23) is its energy-dependent width that is 
parametrized as 

Here To is the decay width when the invariant mass of the daughter combinations is equal 
to rriR. 

The Flatte model is parametrized as 

ArM = — ; j~ ; r. (11) 

m R - s 2 3 - imR{girnP™ + 9kkPkk) 

The constants and gxx are the /o(980) couplings to utt and KK final states respectively. 
The p factors account for the Lorentz-invariant phase space and are given as 




771 (7r" r 7T" 



iil'-^ + Wi-^ (13) 



m 2 (7T" r 7r _ 



For non-resonant processes, the amplitude *4.(si 2 , s 2 3, #,//,/,) is derived from Eq. |4| 
considering that the 7r + 7r~ system is S-wave (i.e. L R = 0, L B = 1) and A fi (s 23 ) is constant 
over the phase space s 12 and s 23 . Thus, it is parametrized as 



P 

A(s 12 ,s 2 3,9ju) = —Jsm 2 e j/ip . (14) 
rriB 
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3.1.2 Detection efficiency 

The detection efficiency is determined from a sample of two million B° — > J/ip{— > 
H + H~)it + it~ simulated events that are generated uniformly in phase space. Both sn 
and S13 are centered at about 18.4 GeV 2 . We model the detection efficiency using the 
symmetric dimensionless Dalitz plot observables 

x = s 12 / GeV 2 - 18.4, and y = s 13 / GeV 2 - 18.4. (15) 

These variables are related to S23 since 

S12 + S13 + s 23 = m 2 B + m 2 j^ + m 2 + + m 2 _ . (16) 



The acceptance in cos 9 j/^ is not uniform, but depends on S23, as shown in Fig. |8j If 
the efficiency was independent of S23, then the curves would have the same shape. On the 
other hand, no clear dependence on s 12 is seen. Thus the efficiency model can be expressed 
as 

£(si2,S23,#//v) =£i{%,y) x £2(^23, 0//vO- ( 17 ) 

To study the cos 9j/^ acceptance, we fit the cos 9 j/^ distributions from simulation in 24 
bins of m 2 (7r + 7r~) with the function 

£2(s ^ 9j ^^ 1 \72a/3 / ^ (18) 



giving 24 values of a as a function of m 2 (7r + n ). The resultant distribution shown in Fig. [9] 
can be described by an exponential function 

a(s 2 3) = exp(ai + a 2 s 23 ), (19) 
with ai = -1.48 ± 0.20 and a 2 = (-1.45 ± 0.33) GeV -2 . 
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Figure 8: Distributions of cos9j/^ for the J/?/>7r + 7r simulated sample in (a) the entire 
dipion mass region and (b) p(770) region. 
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Figure 9: Exponential fit to the acceptance parameter a(s 12) used in Eq. 1 



Equation 18 is normalized with respect to cos 8j/^ . Thus, after integrating over cos 6j/$ , 



Eq. [17] becomes 



£(si2,S23,0j/f)dcos9j/ i j, = ei(x,y). 



(20) 



-1 



This term of the efficiency is parametrized as a symmetric fourth order polynomial function 
given by 

£i(x,y) = 1 + e[(x + y) + e' 2 (x + yf + e' 3 xy + e' 4 (x + yf + e' 5 xy(x + y) 

+40 + vf + ^xy(x + yf + e' 8 xV, (21) 

where the e[ are the fit parameters. 

Figure 10 shows the polynomial function obtained from a fit to the Dalitz-plot dis- 
tributions of simulated events. The projections of the fit are shown in Fig. [IT] and the 
resulting parameters are given in Table [TJ 



Table 1: Efficiency parameters to describe the acceptance on the signal Dalitz-plot. 



Parameter 


Value 




0.142 ±0.010 




0.101 ±0.014 




0.0082 ± 0.0005 




0.027 ±0.007 




0.0052 ± 0.0003 


4 


0.0028 ±0.0010 




0.00074 ± 0.00017 


4 


-0.000105 ± 0.000008 


xViidf 


308/298 
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m 2 (J/\\f ti + ) (GeV 2 ) 

Figure 10: Parametrized detection efficiency as a function of m 2 (7r + 7r~) versus m 2 (<//?/> 7r + ) 
determined from simulation. The z-axis scale is arbitrary. 
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Figure 11: Projections onto (a) i7i 2 (J/iIjtt + ) and (b) m 2 (7r + 7r~) of the simulated Dalitz 
plot used to determine the efficiency parameters. The points represent the simulated event 
distributions and the curves the projections of the polynomial fits. 

3.1.3 Background composition 

Backgrounds from B decays into J/tp final states have already been discussed in Section [2] 
The main background source is combinatorial and its shape can be determined from 
the same-sign 7T tt^ 1 combinations within ±20MeV of the B° mass peak; this region 
also contains the small B~ background. In addition, there is background arising from 
partially reconstructed B° s decays including B Q S — > J/07/(— >■ pj), B° s — > J/i/;(f>(—> n + n^n ), 
and a B° — > J/ipK~7r + reflection, which cannot be present in same-sign combinations. 
We use simulated samples of inclusive B° s decays, and exclusive B° — > J/ip K*° (892) and 
B° — > J/^i^2°(1430) decays to model the additional backgrounds. The background fraction 
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of each source is studied by fitting the J/^7r + 7r candidate invariant mass distributions 
in bins of m 2 (7r + 7c~). The resulting background distribution in the ±20MeV B° signal 



region is shown in Fig. 12 It is fit by histograms from the same-sign combinations and 
two additional simulations, giving a partially reconstructed B° s background of 12.8%, and 
a reflection background that is 5.2% of the total background. 




1 2 3 4 5 

m 2 (jt re) (GeV 2 ) 

Figure 12: The m 2 (7T7r) distribution of background. The (black) histogram with error 
bars shows the same-sign data combinations with additional background from simulation, 
the (blue) points with error bars show the background obtained from the mass fits, the 
(black) dashed line is the partially reconstructed B° s background, and the (red) dotted is 
the misidentified B° — > J/tpK~TT + contribution. 



The background is parametrized as 

B{s 12 ,s 2 3,9ju) = .Bi^^os^) x (l + a cos 2 Oj^) 

where the first part 9 "^ 7r p 7r - converts phase space from su to cos9 nn , and 



(22) 



-Bl(s23,COS 



B 2 (C) 



Pb 
m B 



S2z) 2 + mlTl 



[mt 



x 



1 + Ciq(C) \ cos + c 4 p(() cos 2 6-n 
2[l + Cl g(C)/2 + C4p(C)/3] 



(23) 



The variable ( = 2(s 2 3 — s min ) / (s max — s min ) — 1, where s min and s max give the fit boundaries, 
-82(C) is a fifth-order Chebychev polynomial with parameters hi {i = 1-5), and q(() and 
p(Q are both second-order Chebychev polynomials with parameters q (i=2, 3, 5, 6), and 
Ci, and C4 are free parameters. In order to better approximate the real background in the 
I?g signal region, the J/?/>7r ± 7r =F candidates are kinematically constrained to the B° s mass. 
A fit to the same-sign sample, with additional background from simulation, determines 
6i, q, mo and Tq. Figure [13] shows the mass squared projections from the fit. The fitted 



background parameters are shown in Table [2] 
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The (l + a cos 2 ^j/^ ) term is a function of the J/tp helicity angle. The cos9j^ distri- 
bution of background is shown in Fig. 14 , and is fit with the function 1 + a cos 2 Oj/^ that 
determines the parameter a = —0.38 ± 0.04. We have verified that a is independent of S23. 




Figure 13: Projections of invariant mass squared of (a) m 2 (J/ipiT + ) and (b) m 2 (7r + 7r~) of 
the background Dalitz plot. The points with error bars show the same-sign combinations 
with additional background from simulation. 




cos G 



JA|j 

Figure 14: distribution of the background in cos 6^ resulting from J/ipn+n^ candidate 
mass fits in each bin of cos 9m . The curve represents the fitted function 1 + a cos 2 9ju . 



3.2 Fit fractions 

While a complete description of the decay is given in terms of the fitted amplitudes and 
phases, the knowledge of the contribution of each component can be summarized by 
defining a fit fraction, J 7 ^, as the integration of the squared amplitude of R over the Dalitz 
plot divided by the integration of the entire signal function, 



J a£e<A£(ai 2 , 523,0 



ds\2 ds??. d cos 9 



23 



1 J/i> 



J S (S12, s 23 , ) ds 12 ds 23 d cos 9 



J/0 



(24) 
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Table 2: Parameters for the background model used in Eq. 



23 



Parameter 


T T 1 

Value 


bo 


(4.4 ± 1.2) x 10 6 GeV 


m 


0.767 ± 0.005 GeV 


r 


0.101 ± 0.015 GeV 


bi 


-0.52 ± 0.07 


b 2 


0.22 ± 0.05 


1 


—0.14 ± 0.06 


b{ 


nil _i_ ci ci a 
U.ll ± U.U4 


h 


-0.06 ±0.04 


Cl 


-0.70 ±0.04 


c 2 


-0.4 ±0.3 


c 3 


1.9 ±0.2 


c 4 


0.42 ±0.03 


c 5 


1.7 ±0.8 


c 6 


2.5 ±0.8 


X 2 /ndf 


252/284 



Note that the sum of the fit fractions over all A and R is not necessarily unity due to the 
potential presence of interference between two resonances. If the Dalitz plot has more 
destructive interference than constructive interference, the total fit fraction will be greater 
than one. Interference term fractions are given by 



RR' - -d I / a A afe i( *£ ] A^(s 12 , s 2 3,0j/^)A^' :r (s 12 , s 2 3,0j/4,)ds 12 ds 23 dcos9 



7% = Tie 



and the sum of the two is 



J/0 



fS(s 12 ,s 2 3,Qj/i,) ds 12 ds 23 dcos9 



j/i> 



(25) 



R^R' 

EIE^+E^ R ') = i. (26) 



RR' 



Note that interference terms between different spin- J states vanish, because the d J X0 angular 
functions in are orthogonal. 

The statistical errors of the fit fractions depend on the statistical errors of every fitted 
magnitude and phase, and their correlations. Therefore, to determine the uncertainties 
the covariance matrix and parameter values from the fit are used to generate 500 sample 
parameter sets. For each set, the fit fractions are calculated. The distributions of the 
obtained fit fractions are described by bifurcated Gaussian functions. The widths of 
the Gaussians are taken as the statistical errors on the corresponding parameters. The 
correlations of fitted parameters are also taken into account. 
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4 Final state composition 



4.1 Resonance models 

To study the resonant structures of the decay B° — > J/ip tt + h~ we use those combinations 
with an invariant mass within ±20 MeV of the B° mass peak and apply a J/ip veto. 
The total number of remaining candidates is 8483, of which 3212 ± 80 are attributed to 
background. Possible resonances in the decay B° — > J/ipn + n~ are listed in Table [3J In 
addition, there could be some contribution from non-resonant B° — > J/ipir + 7i~ decays. 

Table 3: Possible resonances in the B° — > J/ip-K+n^ decay mode. 



Resonance 


Spin 


Helicity 


Resonance 
formalism 


/o(500) 








BW 


p(770) 


1 


0,±1 


BW 


o;(782) 


1 


0,±1 


BW 


/o(980) 








Flatte 


/ 2 (1270) 


2 


0,±1 


BW 


/o(1370) 








BW 


p(1450) 


1 


0,±1 


BW 


/o(1500) 








BW 


p(1700) 


1 


0,±1 


BW 


/o(1710) 








BW 



Table 4: Breit-Wigner resonance parameters. 



Resonance 


Mass ( MeV) 


Width ( MeV) 


Source 


/o(500) 


513 ±32 


335 ± 67 


CLEO 


27 


p(770) 


775.49 ± 0.34 


149.1 ±0.8 


PDG 


15 


w(782) 


782.65 ±0.12 


8.49 ±0.08 


PDG 


15 


/ 2 (1270) 


1275.1 ± 1.2 


185.1±1;9 


PDG 


15 


/o(1370) 


1475 ± 6 


113 ±11 


LHCb 


6 


p(1450) 


1465 ± 25 


400 ± 60 


PDG 


15 


/o(1500) 


1505 ± 6 


109 ±7 


PDG 


15 


p(1700) 


1700 ± 20 


250 ± 100 


PDG 


15 


/o(1710) 


1720 ± 6 


135 ±8 


PDG 


15 



The masses and widths of the BW resonances are listed in Table HI When used in the 
fit they are fixed to these values except for the parameters of the /o(500) resonance which 
are constrained by their uncertainties. Besides the mass and width, the Flatte resonance 
shape has two additional parameters g n7T and Qkr-, which are also fixed in the fit to values 
obtained in our previous Dalitz analysis of B° s — > J/ip7i + 7i~ |6|, where a large fraction 
of B° s decays are to J/0 /o (980). The parameters are taken to be m = 939.9 ± 6.3 MeV, 
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g-xK = 199 ± 30 MeV and QkkI 9-k-k = 3.0 ± 0.3. All background and efficiency parameters 
are fixed in the fit. 

To determine the complex amplitudes in a specific model, the data are fitted maximizing 
the unbinned likelihood given as 



N 



125 s 23; ^J/VJ' 



(27) 



where N is the total number of candidates, and F is the total PDF defined in Eq. [T] 
The PDF is constructed from the signal fraction / sig , the efficiency model e(si2, S23, Oj/ip), 
the background model -B(si2, S23, 9j/^), and the signal model S(si2, S23, Oj/^p)- In order 
to ensure proper convergence using the maximum likelihood method, the PDF needs 
to be normalized. This is accomplished by first normalizing the J/ip helicity dependent 
part e(s23,9j/^)Qx(^j/tp) °ver cos 6*^ by analytical integration. This integration results 
in additional factors as a function of S23. We then normalize the mass dependent part 
multiplied by the additional factors using numerical integration over 500x500 bins. 



The fit determines the relative amplitude magnitudes a^ 1 and phases 



defined in 



Eq. [3 we choose to fix a^ 770 ^ to 1. As only relative phases are physically meaningful, one 
phase in each helicity grouping has to be fixed; we choose to fix those of the /o(500) and the 
p(770) (|A| = 1) to 0. In addition, since the final state J/ipn + n^ is a self-charge-conjugate 
mode and as we do not determine the B flavor, the signal function is an average of B° 
and B° decays. If we do not consider tc + 7t~ partial waves of a higher order than D-wave, 
then we can express the differential decay rate derived from Eqs. [3j [4] and [8] in terms of S-, 
P-, and D- waves including helicity and ±1 

dT 



dm^d cos Bawd cos 9 



J/0 

A s s /^ + A s Po e^o cos*™ + A s Do e ir °° (\ cos 2 9^ - \ 



sin 2 Ojfy 



+ 



A s P±i e^ \ sin 9 n7T + A s D±i e l(t>h ^ J \ sin 9 m cos 0„ 



1 + cos 2 9 



j/i> 



(28) 



for B° decays, where At and <f>i are the sum of amplitudes and reference phase for the 
spin-/c resonance group, respectively. The B° function for decays is similar, but 9^+^- and 
Ojfy are changed to n — 9 W + W - and tc — 9j/^ respectively, as a result of using 7r~ and pT to 
define the helicity angles, yielding 

dT 



dm^d cos 9 vn d cos 9 j/^ 

A s So e^o - A Po e^h cos 9 n7T + A^e**^ (\ cos 2 9 m - \ 



sin 2 fl 



J/0 



+ 



A P±l e 9p ^-sm9^-A 



D±i 



i6 s n 

e c ±i 



sin 9 WW cos 9 



1 + cos 2 9 



j/i> 



(29) 
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Summing Eqs. 28 and 29 results in cancellation of the interference involving the A = 
terms for spin-1, and the A = ±1 terms for spin-2, as they appear with opposite signs for 
B° and B° decays. Therefore, we have to fix one phase in spin-1 (A = 0) group (4>p ) and 
one in spin-2 (A = ±1) group (<j) s D±1 )\ the phases of p(770) (A = 0) and / 2 (1270) (A°= ±1) 
are fixed to zero. The other phases in each corresponding group are relative to that of the 
fixed resonance. 



4.2 Fit results 



To find the best model, we proceed by fitting with all the possible resonances and a non- 
resonance (NR) component, then subsequently remove the most insignificant component 
one at a time. We repeat this procedure until each remaining contribution has more than 
3 statistical standard deviation (er) significance. The significance is estimated from the fit 
fraction divided by its statistical uncertainty. The best fit model contains six resonances, 
the / Q (500), / (980), /a (1270), p(770), p(1450), and w(782). 

In order to compare the different models quantitatively an estimate of the goodness 
of fit is calculated from three-dimensional partitions of the one angular and two mass 
squared variables. We use the Poisson likelihood y 2 



28 defined as 



X 



AW 
t=l 



.1 j 



rii + nAn 



ri; 



X j 



(30) 



where rii is the number of events in the three dimensional bin i and Xi is the expected 
number of events in that bin according to the fitted likelihood function. A total of 1021 
bins (iVbin) are used to calculate the x 2 , based on the variables m 2 (J/ipn + ), m 2 (n + 'K~), 
and cos 9 m. The xV n df and the negative of the logarithm of the likelihood, — ln£, of 



7T 1 7T 



the fits are given in Table [5] ndf is equal to iV bin — 1 — iVp ar , where N pai is the number of 
fitting parameters. The difference between the best fit results and fits with one additional 
component is taken as a systematic uncertainty. Figure 15 shows the best fit model 

, m 2 (J/tfj7r + ), cos and m(7r + 7r~). We calculate the fit fraction 
For a P- or D-wave resonance, we report its total fit 
fraction by summing all the helicity components, and the fraction of the helicity A = 
component. The results are listed in Table [6j Systematic uncertainties will be discussed in 
Section pi Two interesting ratios of fit fractions are (0.93loJ2-a23)% for u(782) to p(770), 
and (9.51^ ± 3.0)% for / (980) to / (500). 



projections of m 
of each component using Eq. [24 



The fit fractions of the interference terms are computed using Eq. 25 and listed in 
Table [7j Table [8] shows the resonant phases from the best fit. For the systematic uncertainty 
study, Table [9] shows the fit fractions of components for the best model with one additional 



resonance. 
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Table 5: Values of x 2 / n df and — ln£ of different resonance models. 



Resonance model 


-ln£ 


X 2 /ndf 


Probability (%) 


Best Model 


35292 


1058/1003 


11.1 


Best Model + p(1700) 


35284 


1045/ 999 


15.0 


Best Model + NR 


35284 


1058/1001 


10.3 


Best Model + / (1370) 


35285 


1047/1001 


15.2 


Best Model + / (1500) 


35287 


1049/1001 


14.4 


Best Model + / (1710) 


35289 


1052/1001 


12.6 




Figure 15: Dalitz fit projections of (a) m 2 (7r + 7r _ ), (b) m 2 (J/ip7T + ), (c) cos9j/^ and (d) 
m(ir + Tr~) for the best model. The points with error bars are data, the signal fit is shown 
with a (red) dashed line, the background with a (black) dotted line, and the (blue) solid 
line represents the total. In (a) and (d), the shape variations near the p(770) mass is due 
to p(770) — w(782) interference, and the dip at the K® mass 15 is due to the K® veto. 
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Table 6: Fit fractions and significances of contributing components for the best model, as 
well as the fractions of the helicity A = part. The significance takes into account both 
statistical and systematic uncertainties. 



Components 


Fit fraction (%) 


A = 


fraction S 


ignificance (a) 


p(770) 


cn 0+4.8+2.8 
uz -°-2.9-4.8 


0.63 


+ n 04+ - 06 


11.2 


w(782) 


n p-Q+0. 23+0.27 
u - dy -0.13-0.14 


0.30t°oit ± 0.05 


3.1 


/o(980) 


i cq+0. 77+0.43 
1 - d °-0.50-0.35 




1 


2.5 


/ 2 (1270) 


8.9 ± 1.1 ± 1.0 


0.76 ±0.06 ±0.05 


5.9 


p(1450) 


r q+2.5+5.6 
d -°-1.4-0.9 


0.2 


3+0.17+0.08 
5 -0.13-0.12 


3.2 


/o(500) 


16.2 d= 2.0±|;g 




1 


5.7 


Sum 


95.2 









Table 7: Interference fractions J 7 ^' (%) computed using Eq 
elements are fit fractions defined in Eq [24} 



25 



Note that the diagonal 







p 


u 


P 


fo 


fo 


h 


P 


u 


P 


h 






770 


782 


1450 


980 


500 


1270 


770 


782 


1450 


1270 




\M 




















1 


1 


1 


1 


p(770) 





39.44 


-0.02 


-0.89 























w(782) 







0.18 


-0.05 























p(1450) 









1.47 























/o(980) 











1.53 


2.08 

















/o(500) 













16.15 

















/a(1270) 















6.72 














p(770) 


1 














23.32 


0.29 








w(782) 


1 
















0.41 


-0.07 





p(1450) 


1 


















3.80 





/ 2 (1270) 


1 




















2.14 
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Table 8: Resonant phases from the best fit. 



Components Phase (deg) 



p(770), A = 


(fixed) 


p(770), |A| = 1 


(fixed) 


w(782), A = 


-84 ± 31 


w(782), |A| = 1 


-70 ± 16 


/o(980) 


103 ± 17 


/ 2 (1270), A = 


-87 ± 12 


/ 2 (1270), |A| = 1 


(fixed) 


p(1450), A = 


-162 ±22 


p(1450), |A| = 1 


160 ± 48 


/o(500) 


(fixed) 



Table 9: Fit fractions (%) of contributing components for the best model with adding one 
additional resonance. 





Best 


±p(1700) 


±/o(1370) 


+/o(1500) 


±/o(1710) 


±NR 


p(770) 


62.8l£jj 


59.5+lg 


62.6+^ 


62.4+^ 


63.3+^ 


63.41^ 


w(782) 


o.59±g:g 


0-58lg;S 


o.6o±8:?g 


0.60^ 


0.591^ 


o.59ig:S 


/o(980) 


i -°°-0.50 




i-54±8:H 


i-5o±8:g 


i crcr+0.76 
i -°°-0.51 


i-74l°:l° 


/ 2 (1270) 


8.9 ± 1.1 


8.1 ± 1.2 


8.8 ± 1.1 


8.8 ± 1.1 


8.9 ± 1.1 


8.8 ±1.1 


p(1450) 


5.3lfj 


10.811;! 




4.9111 




4-elri 


/o(500) 


16.2 ±2.0 


15.6 ± 1.9 


16.6 ±2.0 


16.9 ±2.1 


16.3 ±2.1 


21.9 ±3.8 


p(1700) 




3.4±?:S 










/o(1370) 






i q+0.8 








/o(1500) 








± - u — 0.4 






/o(1710) 










0.41^ 




NR 












4.5l? : ? 


Sum 


95.2 


99.6 


96.2 


96.0 


96.7 


105.5 
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4.3 Helicity angle distributions 



We show the helicity angle distributions in the p(770) mass region defined within one full 
width of the p(770) resonance (the width values are given in Table [4]) in Fig. 16 The 
cos 6j/tp and cos 8 nn background subtracted and efficiency corrected distributions for this 



mass region are presented in Fig. 17 The distributions are in good agreement with the 
best fit model. 




300 



1 200 

3 
a 

U 150 



100 



50 



:(b) 




, , | , i i , | , , i , _ 
LHCb : 












, i . , 


- 


1 


-0.5 


0.5 1 

cos 



cos 6 



j/ij. 

Figure 16: Helicity angle distributions of (a) cos9j/^p (x 2 /ndf =15/20) and (b) cos 6^ 
(x 2 /ndf =14/20) in the p(770) mass region defined within one full width of the p(770) 
mass. The points with error bars are data, the signal fit to the best model is shown with 
a (red) dashed line, the background with a (black) dotted line, and the (blue) solid line 
represents the total. 




Figure 17: Background subtracted and efficiency corrected helicity distributions of (a) 
cosOj/^ (x 2 /ndf =20/20) and (b) cos 9^ (rVndf =13/20) in the p(770) mass region 
defined within one full width of the p(770) mass. The points with error bars are data and 
the solid blue lines show the fit to the best model. 
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5 Branching fractions 



Branching fractions are measured by normalizing to the well measured decay mode 
B~ — > J/ip K~ , which has two muons in the final state and has the same triggers as the 
B° — > J/il)n + 'K~ decays. Assuming equal production of charged and neutral B mesons at 
the LHC due to isospin symmetry, the branching fraction is calculated as 

B{B° -> J/^tt+tt-) = ^ s °i^ x B(B~ -> JfyK-), (31) 

Nb- /£b~ 

where N and e denote the yield and total efficiency of the decay of interest. The branching 
fraction B(B~ — > J/tpK~) = (10.18 ± 0.42) x 10~ 4 is determined from an average of recent 
Belle [29] and BaBar [30] measurements that are corrected with respect to the reported 
values, which assume equal production of charged and neutral B mesons at the T(4S), 



using the measured value of p^ go^ = 1-055 ± 0.025 31 . 

Signal efficiencies are derived from simulations including trigger, reconstruction, and 
event selection components. Since the efficiency to detect the J/iIjtt + tt~ final state is not 
uniform across the Dalitz plane, the efficiency is averaged according to the Dalitz model, 
where the best fit model is used. The K® veto efficiency is also taken into account. Small 
corrections are applied to account for differences between the simulation and the data. 
We measure the kaon and pion identification efficiencies with respect to the simulation 
using D* + — > 7i + D°(— > K~7i + ) events selected from data. The efficiencies are measured 
in bins of pi and 77 and the averages are weighted using the signal event distributions 
in the data. Furthermore, to ensure that the p and px distributions of the generated B 
mesons are correct we weight the B~ and -B simulation samples using B~ — > J/ipK~ 
and B° — > J/ipK*° data, respectively. Finally, the simulation samples are weighted with 
the charged tracking efficiency ratio between data and simulation in bins of p and pr 
of the track. The average of the weights is the correction factor. The total correction 
factors are below 1.04 and largely cancel between the signal and normalization channels. 
Multiplying the simulation efficiencies and correction factors gives the total efficiency 
(1.163±0.003±0.017)% for B° ->• J/^vr+vr- and (3.092±0.012±0.038)% for B~ -> J/^K~, 
where the first uncertainty is statistical and the second is systematic. 

Using N B - = 350,727 ± 633 and N so = 5287 ± 112, we measure 

B(B° J/i)n + n-) = (3.97 ± 0.09 ± 0.11 ± 0.16) x 10~ 5 , 

where the first uncertainty is statistical, the second is systematic and the third is due 
to the uncertainty of B[B~ — > J/tpK~). The systematic uncertainties are discussed in 
Section [6] Our measured value is consistent with and more precise than the previous 
BaBar measurement of (4.6 ± 0.7 ± 0.6) x 10" 
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Table 10 shows the branching fractions of resonant modes calculated by multiply- 
ing the fit fraction and the total branching fraction of B° — > J/ip7T + 7r~. Since the 
/o(980) contribution has a significance of less than 3a we quote also an upper limit of 
B (B° ->• J/V>/ (980)) x B (/ (980) ->■ tt+tt-) < 1.1 x 10~ 6 at 90% confidence level (CL); 
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this is the first such limit. The limit is calculated assuming a Gaussian distribution as the 
central value plus 1.28 times the addition in quadrature of the statistical and systematic 
uncertainties. This branching ratio is predicted to be in the range (1 — 3) x 10~ 6 if the 
/o(980) resonance is formed of tetra-quarks, but can be much smaller if the /o(980) is 
a standard quark anti-quark resonance (8l. Our limit is at the lower boundary of the 
tetra-quark prediction, and is consistent with a quark anti-quark resonance with a small 



mixing angle. In Section 7.2, we show that the mixing angle, describing the admixture of 
ss and light quarks, is less than 31° at 90% CL. 

The other branching fractions are consistent with and more precise than the previous 



measurements from BaBar 11 . Using B(lu — > tt + it ) = (1.53_ ' 13 )% |15|, we measure 



B(B° -> J/jju) _ +0 . 24+0 . 31 

T , i nx — U - U1 -0.14-0.16' 

B(B° -> J/ipp ) 

and 

B(B° JAM = (1.5±g3±g5) x io- 5 . 

This is consistent with the LHCb measurement §l§dE^M^ = 0-89 ± 0.19^o!i3, using the 
io — > 7T + 7r~7r° mode 32 



6 Systematic uncertainties 

The contributions to the systematic uncertainties on the branching fractions are listed in 



Table 11 Since the branching fractions are measured with respect to the B~ — > J/-ipK~ 
mode, which has a different number of charged tracks than the decays of interest, a 1% 
systematic uncertainty is assigned due to differences in the tracking performance between 
data and simulation. Another 2% uncertainty is assigned because of the difference between 
two pions and one kaon in the final states, due to decay in flight, multiple scattering, 



Table 10: Branching fractions for each channel. The upper limit at 90% CL is also quoted 
for the /o(980) resonance which has a significance smaller than 3a. The first uncertainty 
is statistical and the second the total systematic. 



Channel 


B(B° J/tpR,R^ 


7T 


+ 71- 


") 


Upper limit of B 
(at 90% CL) 


p(770) 


(0 4O+0.20+0.16N 
1^-^-0.13-0.23/ 


X 


10" 


-5 




w(782) 


(Cy O+0.9 + 1.U 

l z -°-0.5-0.6J 


X 


10' 


-7 




/o(980) 


(6.i±I;^;I) 


X 


io- 


-7 


< 1.1 x 10~ 6 


/ 2 (1270) 


(3.5 ±0.4 ±0.4) 


X 


io- 


-6 




p(1450) 


Icy n +1.0 + 2. 2\ 


X 


io- 


-6 




/o(500) 


(6.4 ±0.8j*j) 


X 


io- 


-6 
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and hadronic interactions. Small uncertainties are introduced if the simulation does not 
have the correct B meson kinematic distributions. We are relatively insensitive to any 
differences in the B meson p and j?t distributions since we are measuring the relative rates. 
By varying the p and pr distributions we see at most a change of 0.5%. There is a 1.0% 
systematic uncertainty assigned for the relative particle identification efficiencies (0.5% 
per particle). These efficiencies have been corrected from those predicted in the simulation 
by using the data from D* + — > tt + D°(-^ K~tt + ). A 0.6% uncertainty is included for the 
J/ijj7r~7r + efficiency, estimated by changing the best model to that including all possible 
resonances. The B° signal yield is changed by 0.5% when the shape of the combinatorial 
background is changed from an exponential to a linear function. The total systematic 
uncertainty is obtained by adding each source of systematic uncertainty in quadrature as 
they are uncorrelated. In addition, the largest source is 4.1% due to the uncertainty of 
B(B~ — > J/tpK~) which is quoted separately. 

The sources of the systematic uncertainties on the results of the Dalitz plot analysis 
are summarized in Table 12 For the uncertainties due to the acceptance or background 
modeling, we repeat the data fit 100 times where the parameters of acceptance or back- 
ground modeling are generated according to the corresponding covariance matrix. We also 
study the acceptance function by changing the minimum IP \ 2 requirement from 9 to 12.5 
on both of the pion candidates. As shown previously (6], this increases the \ 2 °f the fit to 
the angular distributions by one unit. The acceptance function is then applied to the data 
with the original minimum IP x 2 selection of 9, and the likelihood fit is redone and the 
uncertainties are estimated by comparing the results with the best fit model. The larger 
of the two variations is taken as uncertainty due to the acceptance. 

We study the effect of ignoring the experimental mass resolution in the fit by comparing 
fits between different pseudo-experiments with and without the resolution included. As 
the widths of the resonances we consider are much larger than the mass resolution, we 
find that the effects are negligible except for the o;(782) resonance whose fit fraction is 
underestimated by (0.09 ± 0.08)%. Thus, we apply a 0.09% correction to the Co>(782) 



Table 11: Relative systematic uncertainties on branching fractions (%). 



Source 


Uncertainty (%) 


Tracking efficiency 


1.0 


Material and physical effects 


2.0 


Particle identification efficiency 


1.0 


B° p and pr distributions 


0.5 


B~ p and distributions 


0.5 


Dalitz modeling 


0.6 


Background modeling 


0.5 


Sum of above sources 


2.7 


B[B- J/ij)K-) 


4.1 


Total 


4.9 
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Table 12: Absolute systematic uncertainties on the results of the Dalitz analysis. 



Item 


Acceptance 


Background 


Fit model 


Resonance parameters 


Total 


Fit fractions (%) 


p(770) 


±0.9 


+2.0 

— O.l 


+0.6 

— O.Z 


±1.6 


+2.8 

A S 
— 4.5 


w(782) 


±0.08 


+0.23 
—0.06 


+0.11 
—0.10 


+0.028 
— 0.014 


+0.27 
—0.14 


/o(980) 


±0.03 


+0.24 
-0.17 


+0.21 
-0.18 


+0.29 
-0.24 


+0.43 
-0.35 


/ 2 (1270) 


±0.06 


+0.45 
-0.59 


+0.85 
-0.76 


±0.36 


±1.0 


(0(1450) 


±0.10 


+0.5 
-0.6 


+5.6 
-0.7 


+0.4 

-0.3 


+5.6 
-0.9 


/o(500) 


±0.4 


+1.6 
-0.9 


+5.7 
-1.6 


±0.6 


+6.0 
-2.0 






A = fract 


ions (%) 






p(770) 


±1.0 


+ 1.7 

-2.0 


+4.9 
-1.5 


±2.1 


+5.7 
-3.4 


w(782) 


±1.5 


+3.5 
-1.8 


+2.8 
-3.5 


+1.2 

-1.7 


+4.9 
-4.5 


/ 2 (1270) 


±0.3 


±2.4 


+3.7 
-3.4 


±1.5 


±4.5 


p(1450) 


±0.9 


+4.8 
-8.4 


+5.5 
-5.1 


+4.2 
-6.1 


+8.4 
-11.6 


Ratio of fit fractions (%) 


w(782)/p(770) ±0.13 
/o(980)// (500) ±0.3 


+0.41 
-0.11 
+1.5 
-1.1 


+0.18 
-0.16 
+1.0 
-2.1 


+0.034 
-0.022 
+2.2 
-1.8 


+0.47 
-0.23 

±3.0 



fraction and assign an additional ±0.08% in the acceptance systematic uncertainty. The 
results shown in the previous sections already include this correction. 

In the default fit, the signal fraction / s ; g = 0.621 ± 0.009, defined in Eq. [I] is fixed; we 
vary its value within its error to estimate the systematic uncertainty. The change is added 
in quadrature with the background modeling uncertainties. 

The uncertainties due to the fit model include adding each resonance that is listed 
in Table [4] but not used in the best model, changing the default values of Lb in P- and 
D-wave cases, varying the hadron scale r parameters for the B meson and R r esonan ce to 
3.0 GeV -1 for both, replacing the / (500) model by a Zhou and Bugg function 
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and 

using the alternate Gounaris and Sakurai model 35 for p resonances. Then the largest 
variations among those changes are assigned as the systematic uncertainties for modeling 
(see Table 12). 

Finally, we repeat the data fit by varying the mass and width of resonances (see Table ^ 
within their errors one at a time, and add the changes in quadrature. 
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Table 13: Masses of light vector and scalar resonances. All values are taken from 15 
except for the / (500) (27). 



Isospin 


Vector particle 


Vector mass ( MeV) 


Scalar particle 


Scalar mass ( MeV) 







783 


/o(500) 


513 


1 


P 


776 


a 


980 


1/2 


K* 


980 


K 


800 








1020 


fo 


980 



7 Further results and implications 

7.1 Resonant structure 

The largest intermediate state in B° — > J/ipn + -n~ decays is the J/ip p(770) mode. Beside 
the p(770), significant ^2(1270) and /o(500) contributions are also seen. The smaller o;(782) 
and p(1450) resonances have 3.1a and 3.2a significances respectively, including systematic 
uncertainties. The systematic uncertainties reduce the significance of the /o(980) to below 
3a. Replacing the /o(500) by a non-resonant component increases — ln£ by 117, and 
worsens the \ 2 by 192 with the same ndf resulting in a fit confidence level of 1.8 x 10~ 7 . 
Thus the /o(500) state is firmly established in B° — > J/ip7i + n~ decays. 

As discussed in the introduction, a region with only S- and P-waves is preferred for 
measuring sin 2/3 eff . The best fit model demonstrates that the mass region within ±149 MeV 
(one full width) of the p(770) mass contains only (0.72 ± 0.09)% D-wave contribution, 
thus this region can be used for a clean CP measurement. The S-wave in this region is 
(11.9±1.7)%, where the fraction is the sum of individual fit fractions and the interference. 



7.2 Mixing angle between / (980) and / (500) 

The scalar nonet is quite an enigma. The mysteries are summarized in Ref . (36] , and in the 



"Note on scalar mesons" in the PDG 15 . Let us contrast the masses of the lightest vector 



mesons with those of the scalars, listed in Table 13 For the vector particles, the u and p 
masses are nearly degenerate and the masses increase as the s-quark content increases. For 
the scalar particles, however, the mass dependence differs in several ways which requires 
an explanation. Some authors introduce the concept of qqqq states or superpositions of the 
four-quark state with the qq state. In either case, the 1 = /o(500) and the /o(980) are 
thought to be mixtures of the underlying states whose mixing angle has been estimated 
previously (see Ref. [8j and references contained therein). 

The mixing is parameterized by a 2x2 rotation matrix characterized by the angle (p m , 
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giving m our case 



|/o(980)) = cos (p m \ s s) + simp m \nn) 
|/o(500)) = - sin ip m \ss) + cos ip m \nn) , 

where \nn) = -j= (\uu) + \dd)) . (32) 



In this case only the \dd) part of the \nn) wave function contributes (see Fig. |2|). Thus 
we have _ 

_ B(B°^ J/y>/ (980)) $(500) 

Vm b(B° ->• J/V»/ o (500)) $(980)' 1 J 

where the $ terms denote the phase space factors. The phase space in this pseudoscalar to 
vector-pseudoscalar decay is proportional to the cube of the /o three-momentum. Taking 
the average of the momentum dependent phase space over the resonant line shapes results 
in the ratio of phase space factors §g|py being equal to 1.25. 



Using the data shown in Table 10 we determine the ratio of branching fractions for 



both resonances resulting in the ir + n final state as 

B (B° ->■ J/^/ (980)) x B (/o(980) -> tt+tT 
B (B° J/0/ o (5OO)) x B (/ (500) tt+tt 



(9-5±i:I ± 3 



o. 



This value must be corrected for the individual branching fractions of the fo resonances 
into the tt + tt~ final state. 

BaBar has measured B ^j° ^vl7* K + K _ ) = 0.69 ± 0.32 using B — >■ AAA and — > Ktttt 
decays (37). BES obtained relative branching ratios using iJ)(2S) — > 7x c o decays where the 
Xco — > fo (980) fo (980), and either both /o(980) candidates decay into n + Tr~ or one into 

From their results we obtain B ^f e ^1\^ K + ~= ) = 



7T + 7r and the other into A + A pairs 



origin 
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. Averaging the two measurements gives 



ff(/o(980)->A+A-) _ 

B(/ (980)^7T+7r-) - U ' d5 -o-i4- W 

Assuming that the 7T7t and A A decays are dominant we obtain 

B (/o(980) ->• 7T+7T-) = (46 ± 6) %, (35) 

where we have assumed that the only other decays are to 7r°7r°, half of the 7r + 7r~ rate, and 
to neutral kaons, taken equal to charged kaons. We use B (/o(500) — > 7r + 7r~) = |, which 
results from isospin Clebsch-Gordon coefficients, and assuming that the only decays are 
into two pions. Since we have only an upper limit on the J/ip /o(980) final state, we will 
only find an upper limit on the mixing angle, so if any other decay modes of the /q(500) 
(/o(980)) exist, they would make the limit more (less) stringent. Our limit then is 

2 fi(ff>-> J/^/o(980))$(500) 

tan u) m = — 7= r — < 0.35 at 90/o confidence level, 

y B (B° J/0/ o (5OO)) $(980) 
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which translates into a limit 

\<f m \ < 31° at 90% confidence level. 

Various mixing angle measurements have been derived in the literature and summarized 
in Ref. [8j. There are a wide range of values including: (a) using — > it + it + it~ transitions 
which give a range 35° < \(f m \ < 55°, (b) using radiative decays where two solutions were 
found either ip m = 4° ± 3° or 136° ± 6°, (c) using resonance decays from both — > 77r°7r° 
and J/0 —> wkit where a value of <p m ~ 20° was found, (d) using the D ± and Df decays 
into / (980)vr ± and / (980) AT* where values of ip m = 31° ± 5° or 42° ± 7° were found. 

8 Conclusions 

We have studied the resonance structure of B° — > J/07r + 7r~ using a modified Dalitz plot 
analysis where we also include the decay angle of the J/0 meson. The decay distributions 
are formed from a series of final states described by individual 7r + 7r~ interfering decay 
amplitudes. The largest component is the p(770) resonance. The data are best described 
by adding the / 2 (1270), / (500), w(782), p(1450) and / (980) resonances, where the / (980) 
resonance contributes less than 3a significance. The results are listed in Table |6j 

We set an upper limit B (B° ->■ J/0/ o (98O)) x B (/ (980) -»■ n+n-) < 1.1 x 10~ 6 at 
90% confidence level that favors somewhat a quark anti-quark interpretation of the /o(980) 
resonance. We also have firmly established the existence of the J/0 /o (500) intermediate 
resonant state in B° decays, and limit the absolute value of the mixing angle between the 
two lightest scalar states to be less than 31° at 90% confidence level. 

Our six-resonance best fit shows that the mass region within one full width of the 
p(770) contains mostly P-wave, (11.9 ± 1.7)% S-wave, and only (0.72 ± 0.09)% D-wave. 
Thus this region can be used to perform CP violation measurements, as the S- and P-wave 
components can be treated in the same manner as in the analysis of B° s — > J/00 
The measured value of the asymmetry can be compared to that found in other modes such 
as B° — >■ J/ipK° in order to ascertain the possible effects due to penguin amplitudes. 

The measured branching ratio is 

B(B° J/^tt+tt-) = (3.97 ± 0.09 ± 0.11 ± 0.16) x 10~ 5 , 

where the first uncertainty is statistical, the second is systematic and the third is due 
to the uncertainty of B(B~ — > J/ipK~). The largest contribution is the J/0p(77O) mode 
with a branching fraction of (2.49±g;fg±g;||) x 10~ 5 . 
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